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Abstract 

Being motivated by recent experimental studies, we investigate magnetic structures of the Mo 
pyrochlores ^Mx^Or (A= Y, Nd, and Gd) and their impact on the electronic properties. The 
latter are closely related with the behavior of twelve Mo(t2 g ) bands, located near the Fermi level and 
well separated from the rest of the spectrum. We use a mean-field Hartree-Fock approach, which 
combines fine details of the electronic structure for these bands, extracted from the conventional 
calculations in the local-density approximation, the spin-orbit interaction, and the on-site Coulomb 
interactions amongst the Mo(4d) electrons, treated in the most general rotationally invariant form. 
The Coulomb repulsion U plays a very important role in the problem, and the semi-empirical value 
f/~1.5-2.5 eV accounts simultaneously for the metal-insulator (M-I) transition, the ferromagnetic 
(FM) - spin-glass (SG) transition, and for the observed enhancement of the anomalous Hall effect 
(AHE). The M-I transition is mainly controlled by U. The magnetic structure at the metallic 
side is nearly collinear FM, due to the double exchange mechanism. The transition into the 
insulating state is accompanied by the large canting of spin and orbital magnetic moments. The 
sign of exchange interactions in the insulating state is controlled by the Mo- Mo distances. Smaller 
distances favor the antiferromagnetic coupling, which preludes the SG behavior in the frustrated 
pyrochlore lattice. Large AHE is expected in the nearly collinear FM state, near the point of M-I 
transition, and is related with the unquenched orbital magnetization at the Mo sites. We also 
predict large magneto-optical effect in the same FM compounds. 

PACS numbers: 71.70.Gm, 75.25,+z, 71.30.+h, 78.20.Ls 
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I. INTRODUCTION 



The pyrochlores with the chemical formula A2M02O7 (A being the divalent element) ex- 
hibit a number of interesting and not completely understood phenomena,© which present the 
main motivation for the analysis of electronic and magnetic structure of these compounds, 
which we undertake in this work. 

• The pyrochlores A2M02O7 have rather puzzling magnetic phase diagram as a function 
of averaged ionic radius of A-sites, (ta) (Fig. |T])Biii Large (ta) (>-R c ~1-047) stabilizes the 
ferromagnetic (FM) ground state. Typical examples of the FM pyrochlores are Nd2Mo2C>7 
and Gd 2 Mo 2 07. The Curie temperature (T c ) is of the order of 80 K and slowly increases 
with (ta)- Smaller (r a) (<R c ) give rise to the spin-glass (SG) behavior. The characteristic 
transition temperature into the SG state is of the order of 20 K. The canonical example 
of the SG compounds is Y 2 Mo207.i@ Now it is commonly accepted that the exchange cou- 
pling between nearest-neighbor (nn) Mo spins changes the sign around R c and becomes 
antiferromagnetic (AFM) in the SG region. The pyrochlores with AFM interactions in the 
lattice present a typical example of geometrically frustrated systems with infinitely degener- 
ate magnetic ground state.! According to the recent experimental data,i'0 the SG behavior 
itself is caused by the combination of the geometrical frustrations and the disorder of the 
local lattice distortions. The latter is required in order to produce an inhomogeneity in the 
distribution of the nn interactions and freeze the random spin configuration. The origin of 
this behavior itself is very complicated, and we will not be able to address it properly in 
our work. However, our main concern will be in some sense more general: we will try to 
understand which parameter of the crystal structure controls the sign of the nn magnetic 
interactions, and which part of the electronic structure is responsible for the FM and AFM 
interactions in these compounds? Contrary to the widespread point of view that the mag- 
netic ground state is controlled by the Mo-O-Mo angle, we will argue that the key parameter 
is in fact the Mo- Mo distance, which is directly related with the unit cell volume and controls 
the strength of direct interactions between extended Mo(4<f) orbitals. 

According to the experimental data, all Mo pyrochlores which show the SG behavior are 
small-gap insulators. However, an opposite statement is generally incorrect and in different 
compounds the ferromagnetism is known to coexist with the metallic as well as insulating 
behavior J We will show that this trend can be naturally understood in terms of the Mo- 
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Mo distances and the strength of the on-site Coulomb interaction U amongst the Mo(4d) 
electrons. 

• The spin-orbit interaction (SOI) in the pyrochlore lattice should generally lead to a 
non-collinear magnetic ordering associated with the modulation of single-ion anisotropy 
axes.0 The magnitude and the main source of this non-collinearity (whether it is primarily 
associated with the Mo or A sublattices, and whether it reflects the ordering of spin or 
orbital magnetic moments) is largely unresolved problem. In the present work we will try 
to rationalize this question by focusing on the magnetic structure of the Mo sublattice. We 
will show that as long as the system is metallic (that by itself depends on the strength of 
the Coulomb interaction U), the non-collinearity is very small. However, the transition into 
the insulating state is accompanied by the abrupt change of the spin and orbital magnetic 
structures, which become essentially non-collinear. 

• Some ferromagnetic (or nearly ferromagnetic) pyrochlores exhibit rather unusual 
anomalous Hall effect (AHE). A non- vanishing Hall conductivity {cth) emerges near To, 
then monotonously increases with the decrease of the temperature T, and reaches the max- 
imal value (<Tf/~20 fi _1 cm~ 1 for Nd2Mo20y) near T=0. Both magnitude and the temper- 
ature dependence of the AHE is in a drastic contrast with other oxide materials, such as 
colossal-magnetoresistive manganites, for example.0 

The theory of the AHE (originally called as the extraordinary Hall effect) has very long 
and rather controversial historylOB Many controversies were related with the lack of the 
general transport theory and partly - with the lack of good computer facilities in 1950s 
and 1960s, so that many discussions went around very simplified and not always justified 
approximations, such as the free-electron and weak scattering limit, the single-orbital tight- 
binding approach, the classical Boltzmann transport equation, etc.El The modern way to 
address the problem is to use the Kubo formalism and relate au with the real part of the 
antisymmetric off-diagonal element of the conductivity tensor.@ From this point of view 
there are no rigid constraints which would forbid the existence of an, even for periodic 
systems, once the time-reversal symmetry is violated macroscopically.0 

The main question is which mechanism leads to the finite value of a a- The first one was 
proposed by Karplus and Luttinger,El who stressed the importance of SOI in the problem. 
Although the concrete scenario considered by Karplus and Luttinger was not correct, 111 
they were certainly right by arguing that the SOI is essential in order to couple the spin 



4 



polarization of the conduction electrons to the lattice and produce the right-left asymmetry 
of the stationary Bloch states, which may lead to the Hall current, perpendicular to both 
the electric field and the magnetization. 

However, an exceptionally large er# observed in the pyrochlore compounds urged several 
authors to search for a new and unconventional mechanism for the AHE. Such mechanism 
was proposed to be due to the spin chirality (basically, the solid angle subtended by three 
non-coplanar spins) in essentially non-collinear magnetic structure.0El The chirality acts as 
an effective field, which in the combination with the peculiar lattice geometry is believed to 
yield a finite o~h- Although the mechanism itself is rather exotic and there is a great deal of 
cancellations between different elements of the non-collinear spin structure,!!! these authors 
were certainly right by stressing the crucial role of the inter-band transitions in the problem 
of AHE. 

In this work we will go back to a more conventional picture and argue that the large 
value of an in some Mo pyrochlores can be explained by rather unique combination of three 
factors: (i) large SOI associated with the Mo(4oQ states in the nearly collinear FM structure; 
(ii) strong enhancement of the SOI by the on-site Coulomb interactions in the t<i g electron 
systems,^ and (iii) half-metallic electronic structure of the FM pyrochlores.0 

Regarding the Coulomb interaction itself, which in our work will be treated as a pa- 
rameter, we will be able to present a very consistent picture by arguing that C/~1.5-2.5 eV 
in Mo pyrochlores accounts simultaneously for (i) metal-insulator transition; (ii) FM-SG 
transition; (iii) required enhancement for the AHE in the FM pyrochlore compounds. 

The rest of the paper is organized as follows. In Sec. [H]we remind the basic details of the 
crystal and electronic structure of A2M02O7. We consider three characteristic compounds 



corresponding to A— Y, Gd and Nd. In Sec. |TTT| we describe the model Hartree-Fock (HF) 
approach, which combines the ideas of ab-initio electronic structure calculations with the 
physics of degenerate Hubbard model on the pyrochlore lattice. Sees. |V| and [V] summa- 
rize results of calculations without SOI (correspondingly, the basic changes of the electronic 
structure and the orbital ordering, and the behavior of inter-atomic magnetic interactions). 
Sees. |VT| and [VTT] are devoted to the analysis of SOI related properties (the non-collinear mag- 
netic structure, the Hall conductivity, and the magneto-optical effect). Finally, in Sec. |VT l I 
we discuss some subtle points related with the interpretation of experimental data and give 
the summary of the entire work. 
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II. MAIN DETAILS OF CRYSTAL AND ELECTRONIC STRUCTURE 

The pyrochlores A2M02O7 crystallize in a face-centered cubic structure with the space 
group Fm3d, in which A and Mo occupy correspondingly lQd and 16c positions, and form 
interpenetrating sublattices of the corner- sharing tetrahedra. There are two non-equivalent 
types of oxygen sites, and only one internal parameter which may control the properties of 
A2M02O7. That is the coordinate u of the O 48/ sites. 

The single Mo tetrahedron is shown in Fig. [2[ Four Mo sites are located at Ti=(0,0,0), 
T 2=(0, 7, 1), t 3 =(|,0, |), and t 4 =(|, ^,0), in units of cubic lattice parameter a. Each Mo 
site has sixfold O 48/ coordination. The oxygen atoms specify the local coordinate frame 
around each Mo site. Around the site 1, it is given by 

R a/3 = 1 + (1 ~ 8u)6 a0 

1 " V64m 2 -32m + 6' 

where a,/3= x, y, or z. u=j^ corresponds to the perfect octahedral environment. In this 
case Ri=||Rj || is the standard 60° rotation around the cubic [1, 1, 1] axis. u>jq gives rise to 
an additional trigonal contraction of the local coordinate frame. Similar matrices associated 
with the sites 2, 3, and 4 can be obtained by the 180° rotations of Mi around x, y, and z, 
respectively. 

Structural parameters of A2M02O7 are listed in Table |. Corresponding densities of states, 
obtained in the local-spin-density approximation (LSDA) are shown in Figs. |3] and [|. 

In the local coordinate frame, the Mo(4cf) orbitals are split into the triply-degenerate t 2g 
and doubly-degenerate e g states. The splitting is of the order of 4 eV. Twelve t 2g bands are 
located near the Fermi level and well separated from the rest of the spectrum, which consists 
of the broad 0(2p) band spreading from -8.5 to -2.5 eV and either Y(4d) or Gd/Nd(5c?) bands 
located just above the group of the t 2g states. The Mo(e 9 ) states are located in the higher 
part of the spectrum. Interestingly enough that all three compounds are ferromagnetic, even 
in the LSDA, that is rather unusual for the 4d oxides, perhaps except the well know example 
of SrRuO3.il 

The trigonal distortion and different hybridization with the 0(2p) states will further split 
the Mo(t2g) states into the one-dimensional a\ g and two-dimensional e' g representation.!! 
The crystal structure affects the Mo(t2 5 ) bandwidth via two mechanisms (see Table |). 
1. The Mo-O-Mo angle, which decreases in the direction Nd— >-Gd— >Y. Therefore, the in- 
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teractions between Mo(i2g) orbitals which are mediated by the 0(2p) states will also 
decrease. This effect is partly compensated by the decrease of the Mo-0 bondlength for 
A= Gd and Y. 

2. The lattice parameter a and the Mo-Mo distance, which decrease in the direction 
Nd— >Gd— >Y by 2.6%. This will increase direct interactions between extended Mo(4g?) 
orbitals. 

Generally, these two effects act in the opposite directions and partly compensate each 
other. For example, the width of the e' g band is practically the same for all three compounds 
(Fig. (|). On the other hand, the a lg orbitals, whose lobes are the most distant from all 
neighboring oxygen sites are mainly affected by the second mechanism, and the a\ g band- 
width will increase in the direction Nd^Gd— >Y. As we will show in Sec. |V|, this effect plays 
a crucial role in the stabilization of AFM interactions in Y2M02O7, which preludes the SG 
behavior. 

Thus, despite an apparent complexity of the crystal structure, the pyrochlores A2M02O7 
present rather simple example of the electronic structure and in order to understand the 
nature of fascinating electronic and magnetic properties of these compounds, we need to 
concentrate on the behavior of twelve well isolated Mo (£29) bands.S 

III. MODEL HARTREE-FOCK APPROACH 

In this and subsequent sections we will further elaborate this picture, by cutting the 
Mo (£23) bands from the rest of the spectrum, including the on-site Coulomb interactions 
amongst the Mo(4<f) electrons, and solving this problem on the level of HF approximation. 
Hence, our model Hamiltonian has the following form, in the basis of non-relativistic LDA 
eigenfunctions |nk) (LDA is the spin-restricted version of LSDA): 

Wrm'OO = £nk<W + Wf° J (k) + 7# n /0O- ( 2 ) 

The spin indices are already included to n and n', which are running from 1 to 12x2. 
The first term in 7i nn /(k) is the non-magnetic LDA part and e n k are the LDA eigenval- 
ues. 7Y„„/(k) is the matrix of relativistic SOI, which is computed in the conventional way, 
7^^/(k)=(nk||£(r) (L, <x) |n'k), where £( r ) = -r ^ , L is the operator of angular momen- 
tum, cr is the vector of Pauli matrices, V(r) is the self-consistent LDA potential, and c is the 
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velocity of light (c~274 in Ry units). The inclusion of 7i^/(k) into the model Hamiltonian 
is optional, that will be always specified in the text. 

All magnetic effects are caused by the on-site Coulomb interactions, H^ n ,{k), which can 
be written in terms of projections onto an orthonormal set of atomic-like Mo(4d) orbitals 
|Rr/i): 

H% n ,{k) =^^(nk|RT/ii)^ 2 (Rr// 2 |ri / k). 

R.T MlM2 

where R are the translation vectors in the face-centered cubic pyrochlore lattice, and r 
specify different Mo sites in the unit cell. /i={s,m} are the joint indices including the spin 
(s= | or |) and orbital (m= xy, yz, 3z 2 —r 2 , zx, or x 2 —y 2 ) degrees of freedom. Matrix 
elements of the HF potential, U£ lf12 , can be expressed in the most general rotationally- 
invariant form asS 

_ \ (TT _TJ \ (Q\ 

M3M4 

where U^ 1 ^ 3 ^ 2 ^ 4 = (m 1 m 3 \^\m2m 4 )5 SlS2 5 s;iSi are the matrix elements of the Coulomb inter- 
actions. The latter are fully specified by three radial Slater integrals: F°, F 2 , and F A , which 
are assumed to be renormalized from the bare atomic values. Equivalently, one can intro- 
duce parameters of the averaged Coulomb interaction, U=F°, the intra-atomic exchange 
coupling, J=jg(F 2 +F 4: ), and the non-sphericity of the Coulomb interactions between or- 
bitals with the same spin, J B= i i T (9F 2 -5F 4 ). Generally, both U eS =U-J and B contribute 
to the orbital polarization in solids and affect the properties caused by violation of the 
time-inversion symmetry.0 

The density matrix n x =||n^'|| can be obtained from the Green function Q TT i = \\Q^, ||, 
using the standard relation 
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n. 
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where 



Im / deQ TT (e), (4) 



"BZ J ^ — ^nk + 10 



n 



is the spectral representation for the Green function in terms of the eigenfunctions |nk) and 
the eigenvalues e n k of the Hamiltonian (|2|). Q-qz is the volume of the first Brillouin zone. 
Eqs. (§), (|) and (|5|) are solved self-consistently. Then, the total energy is given by 



n BZ t fiifj-2 



lM2 n MlM2 



According to the constraint-LSDA calculations]^ J for Mo is of the order of 0.5 eV, and 
is not sensitive to details of the crystal environment in solids. B can be estimated from J 
using the ratio F 4 /F 2 ~0.63, which holds for the Slater integrals in the atomic limit. This 
yields 5~0.06 eV. The Coulomb U is treated as the parameter in order to consider different 
scenarios, covering both metallic and insulating behavior of A2M02O7. The constraint- 
LSDA calculations for Mo compounds typically yield £/~3.0 eV.0 This value can be further 
reduced by allowing for the (proper) e g electrons to participate in the screening of the on-site 
Coulomb interactions between t 2g electrons.^ 

For the practical calculations along this line, it is convenient to use a short atomic-like 
basis set. For these purposes, e n ^ and |nk) have been first calculated using the nearly- 
orthogonal version of the linearized-muffin-tin orbital method in the atomic-spheres approx- 
imation (ASA-LMTO).il The atomic orbitals |Rr/i) have been identified with the basis 
orbitals of the LMTO method. Atomic spheres radii have been chosen so to reproduce the 
electronic structure obtained by more accurate FLAPW methodJ^EIl especially for t 2g 
bands located near the Fermi level. We have also added a number of empty spheres at the 
8a and 32e positions of the Fd3m lattice. All calculations have been performed on the mesh 
of 12x12x12 k-points in the first Brillouin zone. 

Formally, our approach is similar to the rotationally-invariant LDA+t/,@ except that 
we apply it only to a limited number of bands, directly responsible for the properties of 
A2M02O7 compounds. This will certainly simplify the calculations. In addition, we get rid 
of several artifacts of the conventional LDA+[7 method, such as the substraction of the 
double-counting terms, which is a very ambiguous procedure by itself, and in the number of 
cases leads to a systematic error.il 

IV. CHANGES OF THE ELECTRONIC STRUCTURE CAUSED BY THE 
ON-SITE COULOMB INTERACTIONS 

Let us first discuss results without SOI, assuming FM ordering of the Mo spins. According 
to the LSDA calculations (Fig. f|), the majority (f) spin a lg band is fully occupied and the 
Fermi level crosses the doubly-degenerate e' g band. Therefore, it is clear that at some point 
the Coulomb U will split the e' g band and form an insulating solution with the spontaneously 
broken Fm3d symmetry. Such situation occurs between U= 2.0 and 2.5 eV for all considered 



9 



compounds (Fig. |5], the AFM alignment of the Mo spins may open the band gap even for 
smaller U, that well correlates with the fact that all SG compounds are insulators). 

In the metallic regime realized for small U, the densities of states are similar to those in 
LSDA (Figs. |^ and |j), and we will not show them here again. We only mention that the 
major effect of U is the shift of the |-spin a\ g band to the low-energy part of the spectrum 
relative to the e' g band. Corresponding orbital ordering (the distribution of Mo(4d) electron 
densities) is shown in Fig. || It comes exclusively from the local trigonal distortions of the 
oxygen octahedra and represents the alternating a\ g orbital densities in the background of 
degenerate e' g orbitals. 

Typical densities of states in the insulating phase are shown in Fig. [7|. The a\ g band 
has characteristic for pyrochlores three-peak structure.0 Smaller band gap in the case of 
Y2M02O7 (Fig. H) is related with broad minority ( j) spin a\ g band, whose tail spreads below 
the unoccupied f-spin e' g band. The distribution of the e' g states is very similar for all three 
compounds. Contrary to the metallic state, the orbital ordering is determined not only by 
local trigonal distortions, but also the form of the superexchange (SE) interactions between 
nn Mo sites, and minimizes the energy of these interactions.il Generally, it depends on the 
magnetic state. Two typical examples for the FM and AFM (obtained after the flip of the 
spin moments at the sites 2 and 3) phases are shown in Figs. ^] and respectively. The 
building blocks of the AFM phase are the FM chains of the atoms ... -1-4-1-4-.. and ...-2-3-2- 
3-..., which are coupled antiferromagnetically. As expected for the FM spin ordering ,0 the e' g 
orbitals tends to order " antiferromagnetically" and form two Mo sublattices. Corresponding 
eigenfunctions of the occupied e' g orbitals, obtained after the diagonalization of the density 
matrix in the local coordinate frame, are given by \oi t 4)~a\xy)+(3\yz)+ r y\zx) for the sites 1 
and 4, and (02,3)— a\xy)+ r y\yz)+f3\zx) for the sites 2 and 3, with the coefficients a~0.600, 
/3~0.193, and 7~— 0.769. Clearly, this orbital ordering breaks the Fm3d symmetry: if 
the sites belonging to the same sublattice can still be transformed to each other using the 
symmetry operations of the Fm3d group, the sites belonging to different sublattices - cannot. 
This will generally lead to the anisotropy of electronic properties, including the anisotropy 
of nn magnetic interactions, which will be discussed in the next section. 

The AFM spin ordering coexists with the FM orbital ordering for which a~— 0.774, and 
/5=7~0.442 for all Mo sites, in the local coordinate frame.0 It breaks the Fm3d symmetry 
in the spin sector, but not in the orbitals one. 
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V. NEAREST-NEIGHBOR EXCHANGE INTERACTIONS 



Parameters of exchange interactions between sites r and t' can be calculated asE3 

J TT , = i-Im deTr L [gl^A^gi^e)^} , (6) 
where Gl'^(e) is the block of matrix elements of the Green function with the spin f or |: 

and A x is the magnetic part of the HF potential: 

A T = Ti s {oJ2 t }. 
T?s(L) is the trace running over the spin (orbital) indices. 



Results for nn magnetic interactions for all considered compounds are shown in Fig. [H| 
function of U. We note the following;. 

1. J TT r, which are ferromagnetic for small U, exhibit a sharp drop at the point of transition 
into the insulating state (see Fig. ^| for the behavior of the band gap) . 

2. There is a significant difference between Nd/Gd and Y based compounds: in the Y 
case, the exchange parameters are almost rigidly shifted towards negative values, so 
that the nn coupling become antiferromagnetic in the insulating phase, while it remains 
ferromagnetic in the case of Nd and Gd. 



Similar conclusions can be obtained from the analysis of total energies shown in Fig. |TT|.ti2 
The behavior can be easily understood by considering partial ai g and e' g contributions to 
the nn exchange coupling, calculated after transformation to the local coordinate frame at 
each site of the system (Fig. [TJ). Then, the main interactions are a\ g -a\ g and e' g -e' g . The 



ai g -e' g interaction is small and can be neglected. The interaction between t 2 g and proper e g 
orbitals is of the order of 1.5-1.8 meV and only weakly depends on U and the type of the 
material. Therefore, it can be regarded as the weak FM background. 

Large FM e' g -e' g interaction in the metallic regime is related with the double exchange 
(DE) mechanism, which is the measure of the kinetic energy for the itinerant f-spin e' g elec- 
trons (note that for small U all compounds are either half-metallic or nearly half-metallic - 
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Fig. 0).c2l Some increase of J TT i for small U is caused by the reduction of AFM SE contribu- 
tions, which are proportional to 1/U. As long as the system is metallic, the DE interactions 
are not sensitive to the value of U, and the FM coupling dominates. 

The transition into the insulating state is caused by the localization of the e' g electrons on 
the atomic orbitals. This will reduce the kinetic energy and suppress the DE interactions, 
that explain the sharp drop of J TT i (Fig. [TD| ) . 

The main difference between Y and Nd/Gd based compounds is related with the a\ g - 
aig interactions. Since the |-spin a\ g band is fully occupied and the j-spin band is empty, 
the interactions are antiferromagnetic and the underlying microscopic mechanism is the 
superexchange. Since the SE coupling is proportional to the square of the a\ g bandwidth, 
this interaction is the largest in the case of Y, that explains the AFM character of the total 
coupling realized in this compound for large [7@ 

The anisotropy of the nn interactions, A J= Jyi— Ju, caused by the orbital ordering in the 



FM state is of the order of 1-2 meV (Fig. 10). This is a small value in comparison with the 
sharp drop of J TT i caused by the Coulomb interaction [7. However, it can play an important 
role in the insulating regime, where A J becomes comparable with J TT '. For example, one 
interesting question is whether this A J, of a purely Coulombic origin, is sufficient to explain 
the appearance of the SG state. In addition, AJ has the same order of magnitude as 
the single-ion anisotropy energy (Sec. [VT] ) and the Mo-Gd exchange interaction (Ref. 
Therefore, in the insulating state there is no unique interaction which would control the 
magnetic behavior of the A2M02O7 compounds. All interactions are equally important, that 
represents the main complexity on the way of theoretical description of these materials. 

The AFM spin ordering changes the orbital ordering in the direction which further sta- 
bilizes the nn magnetic interactions,0 and in a number of cases may lead to the charge 
disproportionation within the FM chains.^ For example, [7=2.0 eV appears to be sufficient 
to open the band gap (~0.3 eV) in the AFM phase of Y2M02O7. The solution is accompa- 
nied by the small charge disproportionation, and the values of the spin magnetic moments 
at the four Mo sites are /xf=— /if =1.4 \xb and — /tf=/if=1.5 \lb- Corresponding values of 
the inter-atomic exchange interactions are Ju=— 0.8, Ji3=— 2.9 and Ji4=2.6 meV, i.e. the 
interactions are antiferromagnetic in the bonds 1-2 and 1-3, and ferromagnetic in the bond 
1-4. Therefore, the obtained solution is locally stable. The difference between J12 and J13 
is caused by the charge disproportionation. Further increase of U stabilizes the AFM cou- 
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pling in all Mo-Mo bonds, giving rise to the geometrical frustrations. Typical values of the 
exchange parameters obtained in the HF solution for [7=3.0 eV, which does not show any 
sign of the charge disproportionation, are Ji 2 = Ji 3 =— 10.7 and Ji 2 =— 2.4 meV. 



VI. NON-COLLINEAR MAGNETIC ORDERING 

In this and subsequent sections we will discuss the effects of SOI on the electronic and 
magnetic structure of the Mo pyrochlores. The matrix element of £ (r) calculated in the basis 
of Mo(4g?) orbitals is of the order of £^=120 meV. It defines the characteristic energy gain 
after adding the SOI into the model HF Hamiltonian. The equilibrium magnetic ordering 
depends on the competition between the isotropic exchange interactions, considered in the 
previous section, and the single-ion anisotropies energy. The isotropic exchange interactions 
tend to align Mo spins parallel to each other (at least in the case of the FM coupling). 
However, due to the peculiarities of the pyrochlore structure, there is no unique direction of 
the spin magnetization which would satisfy the single-ion anisotropy energies simultaneously 
for all Mo sites. These anisotropies lead to the distortion of the collinear FM ordering.0 

In the case of FM coupling between nn spins, there are two magnetic structures, which 
are typically discussed in the literature and can be obtained by distorting the collinear FM 
ordering along either (0, 0, 1) or (1, 1, 1) direction (see Fig. [13]). 

1. The "two-in, two-out" structure, where the directions of the magnetic moments 
at the four Mo sites are given by: ei= 7 ^= (r], r], 1— rf), e 2 = 7 ^= (— 77, 77, 1— rj 2 ), 
e 3 = 7 ^= (r),—r), 1—rf), and e 4 = 7 ^= (—r),—r), 1—rf), in terms of the distortion 77 from 
the perfect FM alignments. 

2. The "umbrella" structure: ei=-Wl, 1, 1), e 2 = . 1 (1-277,1+77,1+77), 

' ' y3(l+2ij 2 ) 



Since the Mo(4d) states are strongly polarized by the on-site Coulomb interactions, the 
main contribution to the single-ion anisotropy energy comes from the f-spin electrons. In 
this case there is an one-by-one correspondence between the single-ion anisotropy energy and 
the magnitude of the orbital magnetization, so that the direction which maximizes the mag- 
netization automatically minimizes the magnetic anisotropy energyEl For the t 2ff electrons 
affected by the trigonal distortion, the spontaneous orbital magnetization is determined by 
the following factors:@ 

13 



1. The interaction between occupied a\ g and unoccupied e' g states favors the in-plane ge- 
ometry of the orbital magnetization (i.e., perpendicular to the trigonal axes). The latter 
is proportional to the ratio ^d/Se aig , where Se aig is the distance between the position 
of the dig states and the unoccupied part of the spectrum. Clearly, this interaction will 
vanish when the Coulomb U is large. 

2. On the contrary, the interaction between e' g orbitals favors the easy axes geometry. In 
the pure ionic limit, such solution is the subject for the degenerate perturbation theory. 
Therefore, the wave functions and the orbital magnetization do not depend on £4^. The 
SOI in such situation only lifts the orbital degeneracy and defines the symmetry of the 
e' g orbitals which will be further split by U. The solution is accompanied by the large 
orbital moment of the order of 1 which is parallel to the trigonal axis. 

Thus, when U is sufficiently large, the second mechanism will dominate and one can 
naturally expect the easy axes scenario for the single-ion anisotropics. In the combination 
with the FM coupling between nn spins, it yields the "umbrella" ground state.^1 

The qualitative analysis is supported by results of direct calculations shown in Fig. [14]. 
In the metallic regime, the isotropic double exchange interaction between Mo spins clearly 
dominates over the single-ion anisotropy energies and the canting of spin magnetic moments 
is small (both for the "two-in, two-out" and "umbrella" structure). The orbital moments are 
not particularly strong (<0.1 hb per Mo site). They are canted off the global FM direction 
by about 10°. 

In the insulating state, the isotropic exchange interactions are suppressed and become 
comparable with the single-ion anisotropy energies, which are enhanced. This results in 
the large canting of the magnetic moments, which can be very different for the spin and 
orbital counterparts (in the other words, the spin and orbital moments are not collinear 
with respect to each other) .00 The reason for this difference is the following. Directions 
of the orbital magnetic moments are determined mainly by the single-ion anisotropy terms, 
whereas the interactions between moments located at different sites are typically much 
smaller. Therefore, the orbital moments always tend to align first along the crystallographic 
directions (the trigonal axes, in our case). On the contrary, alignment of the spin magnetic 
moments is strongly affected by the inter-site interactions, which generally lead to a smaller 
canting. 

For large U (~3 eV), the single- ion anisotropies start to dominate over the isotropic 
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exchange interactions. In this limit, the orbital magnetic moments are almost perfectly 
aligned along the trigonal axes (that corresponds to the 55° and 70° canting for the "two-in, 
two-out" and "umbrella" structure, respectively). The value of the local orbital magnetic 
moment is about 0.6 yU#, which is smaller than 1 expected in the pure atomic limit. 

The "umbrella" structure has lower energy (perhaps except the region of very small U). 
The energy difference between the "two-in, two-out" and "umbrella" structures is small in 
the metallic regime (~0.05 meV per formula unit, which corresponds to the coercive field of 
1 Tesla). However, it becomes large in the insulating state, where the interactions between 
the e' g orbitals clearly dominate. 

In the case of AFM coupling between nn spins, the single-ion anisotropies lift, at least 
partially, the high degeneracy of the magnetic ground state in the frustrated pyrochlore 
lattice. The equilibrium magnetic structure realized in model HF calculations for Y2M02O7 
can be obtained by the inversion of the magnetic moments at the sites 2 and 3 of the "two- 



in, two-out" structure (Fig. [15]). For [7=3.0 eV, the values of spin and orbital magnetic 
moments at the Mo sites are 1.44 and 0.51 /zb, respectively. They are canted off the cubic 
(0,0,1) axis by correspondingly 17° and 40°. The AFM interactions are compromised with 
the tetrahedral geometry of the pyrochlore lattice by forming the AFM coupling between 
two of the three projections of the nn spins in each of the Mo- Mo bonds, while the remaining 
projection is coupled ferromagnetically. 



VII. ANOMALOUS HALL CONDUCTIVITY AND MAGNETO-OPTICAL 

EFFECT 

The conductivity tensor, d(u)=d reg (uj)+a D (uj), can be calculated from the eigenvalues, 
£kn, and the eigenfunctions, |kn), of the HF Hamiltonian using the Kubo formula JllEi where 
the inter-band (regular) and intra-band (Drude) contributions are given by (in Ry units) 

/kra' Svn 
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respectively. In these equations, /k n is the zero-temperature Fermi distribution function, 
7r" n ,=(kn|(— iV a )|kn') is the matrix element of momentum operator, u) n i n —e^ n i— £k n , V is 
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the unit cell volume, and r is the phenomenological relaxation time (l/r~68 meV) introduces 
in order to perform summation over the discrete mesh of k-points. The diagonal elements 
of cr D (u;) can be also expressed through the plasma frequencies, Up, as <J ^ a = ^~^~z^rp- 
The conductivity tensor has the form 



a 



^ @xx @xy ^ 

&xy @xx 

\ o zz j 



(9) 



both for the "two- in, two-out" and "umbrella" structure, if the direction of the z axis is 
chosen as (0, 0, 1) and (1, 1, 1), respectively. 

The Hall conductivity, defined as <7H=Reo' X y(0), is shown in Fig. IB as a function of 
U. As expected^ the effect is caused by the inter-band transitions, while all intra-band 
contributions vanish, even in the metallic state, on takes the maximal value around i7=1.5 
eV, which compromises between the value of the orbital magnetization and the metallic 
behavior of the system. On the one hand, U is large enough to unquench the orbital magnetic 
moment of the order of 0.08 hb per Mo site. On the other hand it is small enough in order to 
keep the system in the metallic regime (see Fig. |5[). The factor-two increase of <jh from U =0.5 
to 1.5 eV well correlates with similar increase of the orbital magnetic moment (Fig. |14|) , 
suggesting that both effects have relativistic origin. The behavior coexists with nearly 
collinear FM ordering realized in the metallic regime. The transition into the insulating 
state, which is accompanied by the large canting of the magnetic moments, results in the 
abrupt drop of cr# .@ Therefore, we conclude that the large (jjj observed in the experiment is 
due to the orbital magnetic moments induced by the SOI in the nearly collinear FM stated 
The absolute value of an, ~25 f2 _1 cm _1 , obtained in the calculations is in a good agreement 
with the experimental data for Nd2Mo2O7.0 However, we would like to emphasize again that 
the SOI alone is not sufficient to obtain such a large values of au, and it is essential to 
consider the on-site Coulomb interactions beyond LSDA, which play a very important role 
in the problem and enhance the effect of the SOI.0 

According to our calculations, a xy , in the geometry which specify the form of the tensor 
(§), only weakly depends on the magnetic structure and the direction of the net magneti- 
zation. Large anisotropy of oh observed in the experimentlM! is presumably related with 
the details of the experimental set-up and reflects the " -^-law" . For example, if the electric 
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filed is directed along the (1,0,0) axis, and the current is measures along the (0, 1,0) axis, 
the effect is proportional to a xy /\^3 in the case of the umbrella structure, that roughly cor- 
responds to the experimental situation. Similar arguments apply for the magneto-optical 
Kerr rotation which will be discussed below. 

If our scenario is correct, and large AHE in the Mo pyrochlores can be explained by 
the inter-band transitions in the presence of the SOI (i.e., by the conventional means), it 
is natural to expect similar enhancement for other phenomena related with the SOI. Below 
we consider the complex Kerr effect, which can be also regarded as a possible check for our 
theory in the future. 

The Kerr rotation {Ok) and Kerr ellipticity (e^) can be found from the expression,!! 

tan e K j li6ji 1 - n + 1 + n 



e™« = ^^±^^_ ( 10 ) 
1 + tan 6k 1 + n + 1 — n_ 

which is formally exact and valid for arbitrary Q K and ex- n± = 

\J ^xx — ^— ^£ xy OXG 

the complex refractive indices, in terms of matrix elements of the dielectric tensor, 
E a p(u)=5 af 3+^a a p(uj). 



Results of these calculations for Gd2Mo207 are shown in Fig. |17|, together with the in- 
dependence of the matrix elements of the optical conductivity. We consider the " umbrella" 
structure obtained for [7=1.5 eV, which corresponds to the maximal rotation. Calculations 
for Nd2Mo207 and the "two-in, two-out" structure yield very similar results. We also discuss 
partial contributions of the ai g and e' g states, and show the spectra obtained after removing 
the dig orbitals from the basis.il 

Apart from the HF model itself, we do not use any adjustable parameters. For example, 
the intra-band Drude contribution to the optical conductivity has been obtained in the 
calculations. 

The diagonal part of <t(uS) comes mainly from the inter-band e' g — >e' g transitions. The 
contribution of the a\ g states is rather small. The absorptive (Re) part of <J xx (oS) has a 
maximum around 1 eV. Both the peak position and its intensity roughly agree with the 
experimental data,@ tough the latter reveal an additional structure in the low-cu energy 
part of the spectrum, which can be due to the quantum fluctuation effects,!! beyond the 
mean-field approximation considered in our work.0 

The contribution of the a\ g states to a xy is very important, and can bring up to 50% of the 
spectral weight at certain values of uj. However, the low-u; part of a xy comes mainly from the 
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inter-band e' g — >e' g transitions, which are responsible for almost 90% of the Hall conductivity. 
The large intensity of a xy is due to the combination of three factors: (i) large SOI at the 
Mo sites; (ii) the enhancement of the SOI by the on-site Coulomb interactions, which is 
especially efficient for the degenerate e' g states; (iii) the half-metallic electronic structure of 
Gd2Mo207 and Nd2Mo207: since only f-spin states contribute to the low-a; part of a xy , the 
cancellation between different spin channels does not take place. 

The Kerr rotation spectrum, calculated using Eg. (0), has two structures. The first one is 
close to the plasma edge (ftu;~0.3 eV), and is directly related with the plasma enhancement 
effect.0 The exceptionally large value of 9k (~8°) is also related with large a xy in this 
region, which correlates with the large value of the Hall conductivity. The second structure 
(6k~— 4°) near 0.7 eV is related with the form of R,ea xy (u), which has a deep in this region 
as a joint effect of both e' g —>e' g and ai g —>-e' g inter-band transitions. Smaller U would yield 
somewhat smaller estimates for the rotation angles (correspondingly 4 and —1.5° for [7=0.5 
eV). Nevertheless, the FM pyrochlores A2M02O7 seem to be rather promising materials for 
the magneto-optical applications. 

The anomalous Hall and Kerr effects simultaneously vanish in the AFM phase, due to 
the cancellation of contributions coming from different Mo sublattices in the macroscopic 
sample. Iii 



VIII. DISCUSSIONS AND SUMMARY 



We have presented results of model HF calculations for pyrochlores A2M02O7 (A= Y, 
Gd, and Nd). The model combines details of LDA electronic structure for the Mo(t 2g ) bands 
located near the Fermi level, the spin-orbit interaction, and the on-site Coulomb interactions 
amongst the Mo(4d) electrons. 

The main results can be summarized as follows. 

• The magnetic ground state of the A2M02O7 compounds depends on the Mo-Mo distance, 
which controls the strength of AFM SE interactions between the a\ g orbitals. The latter 
competes with the FM interactions between the e' g orbitals, and dominates when the Mo- 
Mo distance is small. In the frustrated pyrochlore lattice, this corresponds to the transition 
into the SG state, though we cannot specify precisely which mechanism freeze the random 
configuration of the Mo spins. It could be the disorder of the Mo-Mo bondlengths proposed 
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in Refs. ^| and [10]. Another possibility is the anisotropy of nearest-neighbor magnetic in- 
teractions caused by the orbital ordering.il It is also important to note that the Coulomb 
repulsion U should be sufficiently strong in order to open the band gap and suppress the 
FM double exchange interactions between the e' g orbitals. 

• The metal-insulator transition in these systems does not necessary coincide with the 
FM-SG transitions. All metallic compounds are expected to be nearly collinear ferromagnets, 
due to the double exchange mechanism. However, the inter-atomic exchange coupling in the 
insulating state can be both ferromagnetic and antiferromagnetic, depending on the Mo- 
Mo bondlength. The finding is qualitatively consistent with the experimental data.i The 
magnetic structure in the insulating state is expected to be strongly non-collinear, due to 
the competition between the isotropic exchange interactions and the single-ion anisotropies. 

• The large Hall conductivity can be explained by the combination of three factors in 
the nearly collinear FM state: large SOI at the Mo sites, which unquenches the orbital 
magnetization; strong enhancement of this effect by the on-site Coulomb interaction; and 
the half-metallic electronic structure of the FM pyrochlores. We also expect large magneto- 
optical effect in the same compounds. 

A possible way to enhance these properties is to use heavy 5d elements, which have large 
spin-orbit coupling. Unfortunately, all known 5d pyrochlores are non-magnetic, though some 
of them are presumably close to the point of magnetic instability.^ Nevertheless, one can 
try to exploit an old idea known from the physics of intermetallic alloys,il and replace part 
of the Mo sites in A2M0O7 by the 5d elements. Then, the main question is whether the 
remaining Mo sites will be magnetic in the new environment or not. If they will, they will also 
magnetize the 5d states at the neighboring sites, and in this way enhance the non-reciprocal 
characteristics of the considered pyrochlore compounds. As a purely model example, we 
have considered the ordered Gd 2 MoW07 alloy (using the lattice parameters of Gd 2 Mo 2 7 ). 
According to the LSDA calculations, the system is expected to be ferromagnetic. The values 
of (spin only) magnetic moments at the Mo and W sites are 2.1 and 1.5 hb, respectively. The 
situation is highly unusual, because the examples when 4d and 5d elements form a magnetic 
solution are rather rare. Although Gd 2 MoW07 itself may be rather unrealistic example 
from the technological point of view, the idea to use Mo atoms in order to magnetize the 
heavy 5d elements in the pyrochlore lattice seems to be promising. 

In conclusion, the simple HF calculations allow to rationalize many unusual properties of 
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the pyrochlore compounds A2M02O7. It is also important that all these properties can be 
explained by the universal value of £/~1.5-2.5 eV, which is of the same order of magnitude 
as the characteristic bandwidth for the Mo(t2 9 ) electrons. This will certainly place the 
considered compounds to the direct proximity of the metal-insulator transition, meaning 
that the quantum fluctuations beyond the mean-field HF approach will play an important 
role too. This is certainly a very important direction for the extension of the present analysis. 

Finally, we would like to make several comments on the comparison between our calcu- 
lations and the experimental magnetization0111zl and neutron diffractionBIl data. 

(1) The experimental value of the Mo moment (~1.1-1.3 seems to be smaller than 
~1.4 fiB obtained in our mean-field calculations for the nearly collinear FM state, even after 
taking into account the orbital part (~0.1 /i_b), which according to the third Hund's rule 
is antiparallel to the spin one and partly compensate the latter. The small difference is 
presumably related with the quantum fluctuations at the proximity of the metal-insulator 
transition.il 

(2) The experimental data for Nd2Mo2C>7 and SUL2M02O7 typically indicate at the ex- 
istence of two antiferromagnetically coupled sublattices. These sublattices are usually re- 
garded as Mo and A ones. However, according to the electronic structure point of view, the 
spin coupling between these two sublattices should be ferromagnetic,^ being actually consis- 
tent with the experimental data for Gd2Mo2O"7.0 The controversies can be reconciled with 
each other if the local magnetic moment at the Nd/Sm sites has a large orbital counterpart, 
which is larger than the spin one. The AFM orientation of the spin and orbital moments at 
the Nd/Sm sites is again the subject of the third Hund's rule. Then, the magnetization and 
neutron diffraction experiments presumably probe the AFM coupling between spin magnetic 
moments at the Mo sites and the orbital magnetic moments at the Nd/Sm sites. Gd 3+ has 
no orbital moment, and therefore the experimental data for this compound indicate at the 
conventional FM coupling between two spin sublattices. Then, the sharp decrease of the 
magnetization observed in Nd2Mo20"7 at low T,0 which is frequently tempted to connect 
with the anomaly of the Hall conductivity is due to the orbital magnetic structure, not 
the spin one. The change of the orbital magnetic structure seems to be the right subject 
to concentrate at when trying to explain, for example, the temperature dependence of the 
AHE. In this respect, another interesting problem seems to be the behavior of spin and 
orbital magnetic moments at the Mo sites, because already in the ground state these two 
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sublattices from mutually non-collinear magnetic structures. Therefore, it is reasonable to 
expect very unconventional spin and orbital dynamic at the elevated temperatures. 
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TABLE I: Structural parameters of A2M02O7 (A= Y, Nd, and Gd): the cubic lattice parameter a 
(in A), positions of 48/ sites [it, |, |] (units of a), the distances Mo-Mo and Mo-0 (in A), and 
the angles Mo-O-Mo and O-Mo-0 (in degrees). 



compound 


a 


it 


^Mo-Mo 


dtvIo-O 


ZM0-O-M0 


ZO-Mo-0 


Y 2 Mo 2 7 a 


10.21 


0.33821 


3.6098 


2.0171 


127.0 


99.5 


Gd 2 Mo 2 (V 


10.3356(1) 


0.33158 


3.6542 


2.0123 


130.4 


97.2 


Nd 2 Mo 2 7 c 


10.4836(2) 


0.32977 


3.7065 


2.0332 


131.4 


96.6 



a Ref. I a is from Ref. [5] (corresponding to 4K). 
6 Ref. § 
c Ref. |. 
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FIG. 1: Experimental phase diagram of j4 2 Mo 2 Oy: transition temperature against the averaged 
ionic radius. The values for Nd 2 Mo 2 07 are taken from Refs. || ((g)) and [2^ (0), for Sm 2 Mo 2 07 
- from Refs. | (+) and || (x), for (Smi_ x .D y:r ) 2 Mo 2 07 (x= and 0.2), (Ndi-^Tb^Mo^ 
(x= 0.5), (Smi_ 3; Tb ;c ) 2 Mo 2 07 (x= 0.5), (Gdi^Dy^Mo^ (x= 0, 0.1, 0.2, 0.4, and 1), and 
(Gdi_ a; Tb ;c ) 2 Mo 2 07 (x= 0.25, 0.5, 0.75, and 1) - from Ref |, and for Y 2 Mo 2 7 - from Ref. [7|. 
Dotted line shows the boundary between normal ferromagnetic (right) and spin-glass (left) behav- 



ior. 
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FIG. 2: Relative position of the Mo and O 48/ sites in pyrochlore lattice. 
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FIG. 3: Total and partial densities of states of Y2M02O7 and Gd2Mo207 in the local-spin-density 
approximation. The Mo(t2 9 ) states are located near the Fermi level (chosen as zero of energy). 
The Mo(e 9 ) states emerge around 4 eV. 
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FIG. 4: Mo(t2 g ) states in the local coordinate frame, split into the one-dimensional a\ g and two- 
dimensional e' g (denoted as e g (t2 g )) representations by the local trigonal distortion. The Fermi 
level is at zero. The contribution of proper e g states to this region is small. 
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FIG. 5: The band gap as a function of Coulomb U. Here, the band gap is defined as the distance 
between the bottom of 9th and the top 8th band, so that for C7<2.0 eV there is an overlap between 
these two bands, corresponding to the negative value of the band gap, while U>2.5 eV opens the 
real gap. 
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FIG. 6: Orbital ordering in the ferromagnetic phase of Y2M02O7 obtained for C/=0.5 eV. 
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FIG. 8: Orbital ordering in the ferromagnetic phase of Y2M02O7 obtained for [7=3.0 eV. Two 
orbital sublattices are shown by black and grey colors. 
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FIG. 9: Orbital ordering in the antiferromagnetic phase of Y2M02O7 obtained for [7=3.0 eV. Two 
spin sublattices are shown by black and grey colors. 
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FIG. 10: Nearest-neighbor exchange interactions calculated in the ferromagnetic state. The orbital 
ordering realized for large U destroys the equivalence of some of the Mo-Mo bonds, that generally 
leads to the inequality J\2=Jiz^ z Jia- Two such parameters, J±2 and J14, are shown by dark and 
light symbols, respectively. 
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FIG. 11: Total energy of the antiferromagnetic state relative to the ferromagnetic state. 
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FIG. 12: Contributions of a\ g and e' g orbitals (coming from the t2 9 manifold in the local coordinate 
frame) to the exchange constant J14. 
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FIG. 13: Two magnetic structures realized in the case of the ferromagnetic coupling between 
nearest-neighbor Mo spins: "two-in, two-out" (left) and "umbrella" (right). The relative orienta- 
tion and the magnitude of the spin (light vectors) and orbital (dark vectors) magnetic moments 
corresponds to Gd2Mo2 0y for two different values of the Coulomb interaction U. 
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FIG. 14: The angles formed by the spin and orbital magnetic moments and either (0, 0, 1) ("two-in, 
two-out") or (1,1,1) ("umbrella") axes; the absolute value of the orbital moment; and the total 
energy of the "umbrella" structure relative to the "two-in, two-out" structure as a function of 
Coulomb U. Filled and open symbols correspond to the "two-in, two-out" and "umbrella" struc- 
ture, respectively. In the case of "umbrella" structure, there are two different magnetic moments: 
the ones at the site 1 are shown by regular symbols, O an d and the ones at the site 2, 3 and 4 
(forming the top of the "umbrella") are shown by crossed symbols. The spin magnetic moment is 
of the order of 1.5-1.6^, and weakly depends on U. 
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FIG. 15: Magnetic structure in the case of antiferromagnetic coupling between nearest-neighbor 
Mo spins. The relative orientation and the magnitude of the spin (light vectors) and orbital (dark 
vectors) magnetic moments corresponds to Y2M02O7 for U =3.0 eV. 
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FIG. 16: Hall Conductivity in Gd2Mo207 and Nd2Mo207 as a function of Coulomb U . Filled and 
open symbols correspond to the "two-in, two-out" and "umbrella" structure, respectively. 
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FIG. 17: Matrix elements of optical conductivity and the Complex Kerr effect for Gd2Mo2C-7 
corresponding to the "umbrella" ground state realized for U =1.5 eV. Black lines show total spectra. 
Gray lines are obtained after excluding the a\ g orbitals from the wave functions (and roughly 
correspond to the partial contributions of the e' g bands). 
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